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Tlus work provides analytical survey pf current concepts regarding 
electrostatic ion-cyclotron turbulence (theory and experiment) , and re- 
garding inclined double potential layers in the magnetospheric piasma. 

It examines the questiph of anomalous resistance governed by electro- 
static ion-cyclotron turbul,ence, and one-dimensional and two-dimensional 
models of double electrostatic layers in the inagnetospheric plasma. 
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ANNOTATION 


This work provides an anaiyti^i:al survey of current concepts regarding electro- 
static ion-cyclotron turiaulence (tl'isory aid ej^riment) , and regarding inclined 
double potential layers in the magrv^?tos£iieric plasma. 

It examines the questicxx of ancwfaious resistance governed by electrostatic ion- 
cyclotron turiiulence, and one-dimensioiial and tvro-dimensional models of double elec- 
trostatic layers in the magnetospheric plasma. 
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lON-CYCLOTRCN TURBULENCE AND OBLIQUE 
DOUBLE LAYER IN MAGNETOSPLiERIC PLASMA 

Ey 

V. A. Liperovskiy, M. I. Pudovkin, G. A. 

Skuridin, and S. L. Shalimov* 

This survey reflects the rising interest in studying electrostatic ion-cyclo- /3** 
tron (EIC) turbulence vhich developed because of a number of new measiuretnents made 
in the near-earth space plasma on the S3 - 3 satellite. 

In addition to experimental proof of the ElC-wavas in the magnetospheric plasma, 
the survey examines the main conclusions of the linear theory of ElC-instability. 

A detailed stud^ of it is found in the pioneer works of Drutraiond and Rpsenbluth [1] , 
and Lominadze and Stepanov [2] . 

Study of the nonlinear stage of ElC-instability evolution was made by Petviashvili 
[3], Dum and Di^ree [4] and others. Sagdeyev [6, 7, 8], Sagdeyev and Galeyev [9], 
as well as the book of Artsimovich and Sagdeyev [10] eiramined the prcAilati of arsomf- 
lous resistance governed by plasma turbulence, Kindel [11] , lonson [12] , lonson ^ 
et al. [13] Dakin et al. [14] and Hudson et al. [15] examined a similar problem 
in relation to EIC-turiDulence. 


*Institute of ^saoe Research of the USSR Acadeny of Sciences . 
**Numbdrs in margin indicaije pagination in original foreign t«^. 
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otie appropriate experimental laboratory studies of ElC-turbulence were 
made. Ihe most important has become the experiments of Dakin et al. [14] , 
Correll et al. [16,17] , and Bohmer and Pornaca [18] . As applied to longi- 
tudinal currents in the auroral zcme of the magnetosphere, the role of the 
EIC - instability was first established Kindel and Kennel [19] , /4 

found that, among the basic three wave instabilities excited in a plasma 
with a current (Buneman, ion-sound, and EIC), the ElC-instabi] ty has a 
lower excitation threshold in the magnetosphere. 

In 1963 Petviashvili showed that within the framevork of one-dimensional 
quasilinear theory , in an infinite, uniform, magnetoactive plasma the 
occurrence of a plateau on the electron distribution function stabilizes 
the ElC-waves [3] . The effective frequency of electron collision with 
waves thus equals 

where ■Ou are ion-electron collisions; u - drift velocity of electrons 
along the magnetic field which, for EIC-i.nstability, has a much smaller 
thermal v^. Such a small effect of collective collisions must naturally 
be disreg^ded. . 

Thus for a long time it was assumed that ElC^'instability, in spite of 
a lower excitation threshold as compared with the ion-round one and in 
spite of the absence of the requirement for non-isothermality, was not 
important in the problem of anomalous resistance in the magnetosphere. 

It was confirmed by other authors - Vedenov [20] , Lysak [22] , Dum and 
Dupree [4] that collisions, three-dimensional effects and the presence of 
longitudinal constants of the electric fields can overcome the require- 
ment for a"plateau", and saturation of the turbulence level is performed 
at a higher level (for example, under the influence of the mechanism of 
expansion of ion resonances [4]). As a result, the anomalous resistance 
caused l:y the ElC-turbulence increases. Experiments on the satellite S 3-3 
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[15] also showed that a "plateau" is not formed. Later in apace experi- /5 
ments it vms found that the level o| ElC-turbulence, connected with 
longitudinal currents was so great ^iihat the frecpier^y of collective 
collisicxw of particles with waves leeds to a large rjioroalous resistance 
in the longitudined. currents. Ihus the study of the 2 unomalous resistance # 
caused by EIC tiurbuleiy^e, is in|»rtant for understanding nany large-scale 
geophysical processes. 

•flie study [i\0] presents data on laboratory e}^)eriments indicating that 
the develc^ment of strong wave turbulence of the ion-sound type leads to 
the occurrence of direct potential layers in the current-carrying plasma, 
i.e., layers in vAiich the field is parallel to the magnetic field. Ohe 
double layers, which develop in the plasma at a very high ElC-turbulence 
level, are inclined with respect to the magnetic field (so-called oblique 
layers of tile IC-type) . 'flie physical model of this inclined layer was 
developed in a series of articles fcy Swift [51-55] . 

The basic features of this model are discussed below, 

1. EIC-INSTABILITO; EXPERIMENT AND THEORY 

1.1. EXPERIMEfiTHi DISCUSSION OF EIC-WAVES IN THE 
MAGNETOSPHERE 

Hiere have been very few d^vect observations of ElC-waves in the 
magnetosphere up until recently. 

EIC waves were observed in the region of diurnal polar cusps oA the 
satellite OGO-5 [42] during a strong magnetic storm on 1 November 1968. 

A frequency was observed which corresponded to the intensity maximum of 
waves located between 0.67 and 0.87 

In a rocket experiment [25] at an altitude of about 600 km above the 
Earth close to the au^ral arc boundary strong electrostatic fluctuations 
were determined about (at frequencies of 40 Hz) . The 


a 


# 


* V 


sijnj^taneously obsecved strong heating of the lc»m is associated the /6 
authbirs of [25] with the presence of EIC instability in these regiois. 

Further, the electrostatic wave turbulence was studied in detail in 
[34,45], vhere a fairly high le(!/el of wave energy was noted at frequencies 
close to Ihese ejcperimenti,^, however, did not nake a special study of 
the EIC waves. American works in the last several years from 1975-1978 
which cover EIC turbulence in the cosmic plasma are studied in [43] . 

An important step forward towards stud/ing wave tiurbulence in the mag- 
netosphere was the studies of the Mozer group which were made based on 
measurements of the S 3-3 satellite [15,24], Obese pubi.icaticns give the 
results of measurements (at altitudes about in the aviroral zone) of 
intensity of both the constant and variable 6n of the plasma. Obe measure- 
ments showed the presence of regions of strong longitudinal electrical 
fields and wave electrostatic turbulence. Obese regions were recorded 
at varying local time (M) on L-shells from L'=4.7toL>6. Oypical 
dimensions of the regions were about 35 kilometers.. Obis corresponds 
to 10 km in the ionosphere. 


Figure 1.1 [24] presents an example of the measurements of spectral 
powers (6n/n)^ and (6E)^ for the 217th day of 1976. It should be 
noted that, in addition to these main measurements, measurements were 
made (evidently with much less accuracy) of the magnetic field. 

Obey made it possible to establish a rough estimate of the magnitude of 
the density of the longitudinal currents. 
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average or typical values of tlie measured paramsters are the followings 

^9 ^ » *5 I 'V/H^ 

j" ■ ' TTiOO ^ I 

, ft 

Analysis of the findings pemitted the authors of [24] to assert that the 
obs^^rved turbulence is ion-cyclotron which correspcaids to the Drutmiand-Rosenblutli 
mode. iMs conclusion was drawn based on the following arguments. First of all, 
in the spectra presented in Figure 1,1, tliere are clear nvaximums of spectral power 
at frequencies slightly greater than oigi, 2wg£, 3wg^ as should be according to the 
linear tlieory of ElC-instability, 

In one of the typical exanples « 1.4, and this , according to the liiiear 
theory (see Figure 1.4) corresponds to the ratio Te/Ti scmewhat greater than one^ 
Thus, altliough direct measurements of Te and Ti were not ^lade, certain information 

regarding the ratio of tenperatures Wrsr? extracted by the authors of [24] from tile 

\ 

set of experimental data and the hypothesis regarding the correctness of tlie linear 
tlieoty of ElC-instability for this case. Hcwever, estiiiates for tenperatures can 
hardly be considered ejperimental fact, since currently it is far from clear how 
much the conclusions of linear theory can be used to analyze tiie turbulent quasi- 
stationary state. 


Secondly, in certain cases (vhen the time for passage of the region of wave 
electrostatic turbulence was greater than the period of rotaticm of the satellite 
around its own axis) , polarization of Icw-frequency oscillations was defined. It 
was found that intensity of the variable electrical field is the maximum in the 
direction perpendicular to the main magnetic field Bq with accuracy t 15® . This 
conclusion is an inportant confirmation of tlie fact that the observed waves corres- 
pond to the ElC-waves. 
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'me authors of [24] estimated tlie tai|?erature of the electrons as follows. 
Since tlie quantitieB of fluctuations in tlie electrical field 6 e and density of 
electrons $n for slow waves are associated the ratio 

tKa»ffafiwT» t»wT /9 

(wheira k « wave number) , then by using measurements of ki and iSE one can determine 
tlte product of the wave nuntoer k times TgCkTe" 5 x 10 “^ eV/m) . Hie autliors of [24 ] , 
taking into oo.^sideratian the Dc^pler effect in measuring frequencies frcnv a moving 
satelliter based on the results of a tlieory developed in [46] , estimated the jrfiase 
velocity of waves and their wave number (k~ 1.6 x lO*'^ m'^^) . This permitted fur- 
ther estimate of the tenperature (Tq« 3.4 eV) . The conclusion that the observed 
wave turbulence corresponds to tlie Dnarmond^^Rosenbluth mode is also confirmed by 
the fact tiiat with direct si^titutiQn of the ofetaiiied estimates for wave parameters 
in the actual part of the oorresponding dispersion equation (see equation 2.1 of 
this survey) it is satisifed with a fair degree of accuracy. 

It should be noted that these studies did not iraJ;i measurements of fluctu- 
aticms in the magnetic field/ i'.e . , it was not experimentally proved that the waves 
only have electrostatic naturel 

A shortocmiiig of the discussed reservations is the fact that measurements of 
the variable electrical field were made only to frequencies 5,500 Kz (and in [44] 
data are presented to frequencies 5 1000 Hz ) , while the plasma ion frequency Wp^s 2 
X lO^ Hz could be unnoticed because of the ion-^sonic turbulence vMch possibly 
exj,sts simultaneously with the cyclotron. 

Hius/ the set of wave experiments made on the satellite S 3*53 yields rich 
experimental information necessary to analyze the nature of tlie electrostatic 
wave turbulence in tlie magneto^here. Ihe cycle of works by the Mozer grovp to 
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investigate phenoiiii?m of anonalous resistivity , double layers and electrostatio 
shockwaves baaed on tliese exjex^Jnents and a nuntoer of ^.y^th©ses draws tie conclu- 
sion that the E^C-instability plays a definite role in these giienciitena. However, 
the set of ejqperimfiBntal data, frcm our viewpoint, is still not sufficient for 
reliable confirmation of the definitive role of the ElC-iiode in Idle problem of 
anomalous resistivity in the raagnetosEdiere. One can only assert that the EIC- 
turbulence really exists. 


1.2. MAIN C30NCUJSI0NS OP TOE IiTOEAR 
TOEORK OP EIC-INSTABIliITO 


If one considers that the tenperature of ions and electrons Tg and for the 

plasma in tlie magnetic field has a finite value vAiich differs from zero, then one 

can find the ifellesdng effect [1, 2| . It is found that with 0 sn/2 { 0—angle between 

direction of external magnetic field % and wave vector k) a number of new branches 

of natural plasma oscillations cppear whose frequencies with k*^ 0 and k**-" a£proadi 

' «» *• 

the frequency n[ojBo | < (A* | h . TOese oscillations are 

called cyclotron. 

TOe behavior of the cyclotron oscillations significa.itly depends on the proximity 
of tlie angle 0 tO’Tr/2. TOerefore cases of transverse 

(o* I , Itimxoe\0«o, » TiVk,, , '»f* • 

and quasitransverse dissemination of cyclotron oscillations are isolated. TOe condi- 
tion of quasitransverse nature is written asj « u) V*ICn . 


Vte now assime that low-pressure plasma is examined for which tlie following oondi- 
tions are fulfilled (^plied, in particular, to the magnetosphere) 

« COS^d « 1. 



It new turns out (see for exanpl© (271 ) , that front the lon-cycXotron oscillations 
it is possible to isolate almost electrostatic (potential, ’Xongituclinal) ion- 
Oyclotron oscillations. If the latter oonditic«s are not fulfilled, then the 
oscillations are not potential. In order to obtain electrostatic ion-<^lotron /II 
oscillations in the case of their quasitrans verse dissemination, we will exanvine 
the dispersion equation for the Maxwellian functions of distribution which des- 
cribe the longitudinal plasma oscillations in the magnetic field [27] 5 

£• m HI t Cc * 0 

^ -tf** « yiTT^ I 


Ij^(Ma) — modified Besseiian function. The function W is called tl>e integral of proba- 
bility from the conplex argunent or Grant's function. In tlie case viiere the argu- 
ment of the W function is great, or small as cc5npared to a unit, this function can 
bfi approximated by the following approximate formula which follow from its deter- 
muvs^cion [36] : 

•> 

f 4 ♦ % 

V </ia »» 

tw) |r| « 4 

^0) • ^ ^ “ 

miy 
\] 

Only the real arguments of the W function are further examined everyvhere, NW'7 in 
order toisolate the longitudinal ion-cyclotron oscillations in the case of their 




quasitrajTsverse dissemination, we add the following to the hypotheses and conclu- 


sions made: 


iCi . 
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Tlien, after using the asyirptotic facmulas for W(x ) , tlie contributions of and 
in the dispersion equation for iongitaidinal osciUaticsis ^ 6^*0 will 

be determined by the fontulas; 

«»• (1-.2) 

r,(^n) > (.(*■) , f* • “hi* 


Only one antiHermitian term responsible for the Cerenkov absorption of oscillations 
by the electrons is left in the expression for since according to the condition 
*»> Kh’l* , the electron cyclotron dampening of the examined oscillations 
has an exponential order of smallness. The condition means tliat 

Kh'di /0%i i< ** i.e., dissemination of the wave occurs almost perpendicularly to 
the magnetic field. At the same time, weak ion cyclotron damping of the oscillations 
is guaranteed,. We will now assume for simplicity that I ' ' '* • , when one 

can ignore in the dispersion equation the unit as ocnpared to and Sg. Then the 

equation for tl^ real part of the c^cillation frequency v*))f 

in our case adopts the following appearance; 



From here we obtain: 


4 4 


z 

n * 0 


fSllLl 4 


' M ) 


( 1 . 3 ) 
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Sincje ' , then ^ ‘ J 

It is cxsnvenient to solve the equation ' graphically 

♦ 

(see Figure 1.2 [27] ) . The solutions to *0 * H > It 8 *»* jcorresponds to (he 

points of intersection of str^^t line y = 1 + /T^ am' the curve y - f (to) . 



Figure 1.2. Graphic Resolution of the Disper- 
sion JSquation 


/M 


It is ^parent frc«n Figure 1.2, that with a rise ixi V^/Tq, lo^"'^ (k) approaches 
Hie behavior of the frequency (k) depending on "k” is sdiematically 
shown in Figure 1.3 [27] . It follows frocn Figure 1.3, that for long (kpj^« 1) 
and short (kp^ ^"^l) waves, the frequencies w(m) (k) close to nwg^^. Iherefore one 


can assume ^ - 

Ivdiere h order to find ^Ck)t we return to equation (1.3). Ey 

substituting instead of “ its ejqpression (1.4) ,‘ we bring, (1.3) to the follcwing 


agpearance , 4 ♦ H 

7. n 
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Figure 1.3. Schematic Dependence of Frequencies 
of Longitudinal Ion-Cyclotron Oscil' 
lations on Qiantity of Wave Vector 


Ercm here we find . 

( 1 . 5 ) 

Vfe will lindt ourselves in tlie further presentation only to the first harmonics 
(n = 1) . This limitetion is justified ty the fact tliat the weakly-danpening oscil- 
lations with a large nurrber of the branch ”n" do not exist [27, 16] . 


Ihen for the first liarmonifjs instead of (1.4) we have 


( 1 . 6 ) 


vhere 


■ j- I 


r.{r) 

^ ^ ■ Ik 4 "At Cr 


If we now total the infinite sum for "n," ty using the properties of the functions 
r.(r) : " and the recurrent ratios for the Besselian function, 

then 




i + Tt/l* - 


(1.7) /16 




vdiexe ^Cf*) a ji**"^^** » O ®Cf^) < ^ . We will examine the area 

of aH>licabilitY of formulas (1.6) and (1.7). 

In the case of T£» Tg, i.e., in plasma with hot ions and cold electrons (wiiich, 
for exairple, occurs in the transitional region of the magnetosphere and in the area 
of diurnal polar cusps) we have 

t ^ if)) • (1.8) 

Vfe note that this formula expresses the proximity of (k) to but now already 
witli randan values f* f '‘.fj*)* (^)) . Derivation of this usually employed 

formula is not equivalent to the result that is obtained if it is assumed y = 0 (then 
G(o) =1, 1 - G = 0, but T 2 (o) = 0 and u = Wgj|^) . in the case of T^/T s i for the 
maximum quantity = 0.2 wiiicli is reached with M= l.S [19], it is easy to see that 
Aj <0.2, therefore w = 1.2(^£. When Tg/Tf »1 ard 0 (so that G(M )"*• 1) , the cor- 
rection Aj^ can beocme large, since under these conditions / ^ ^ • * 

Ihe numerical solution to equation (1.3) shewed tliat the relationships for different 
ratios of tenperatures have an sppearanoe presented in Figure 1.4 [27] . Ihe elec- 
trostatic ion-cyclotron fluctuations oscillate, for exanple, with drift of tiie elec- 
trons in relation to the ions. If the plasma electrons vmidi have Maxwellian dis- 
tribution for velocities drift under the influence pf an (external electric field 
along tlje external magnetic field with velocity u, thcjn in the electron term 

in the dispersion equation (1.2) w should be replaced by w - k^u. The f requency 
of tlie longitudinal oscillation with u is determined ty formula (1.6) as 
before. 


The ejq)ression for the increment in the oscillations in the linear theory is /18 


found from the formula, 


y Tim €• ( ff ) 

W* • , Xu 1 IiM toC«) 




Figure 1.4. Dependency of Frequency and Longi- 
tudinal lon-cydotron Oscillations 
on (k2p^2) for Different Values of 
the Parameter Te/Tf. Tlie values 
Tg/Ti ec^l to 0.1, 0.3, 1 and 3 
respectively correspcaid to curves 
1-4. 


Hierefore, taking into consideration the small antiHermitian terms in (1.2) we 
obtain for the first harmonic 


(1.9) 


Vfe new find from (1.9) an estimate for the critical drift velcjcity xig vMcdi deter- 
mines the thre^old for development of the instability. In this case, y= 0, there- 
/ ^ U.6), and by 

Ij designating S* » we find that /6i , tlierefore 
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By minimizing the function K a) in relation to k„, kjj^, i.e., we solve the 

equations . n/ As a result, we fi^g^at Ifif. ^ ^ (so that 

Vi = 0.2) and (uJ » 2 f ’ 

With these values ti/it* s-fO and 

. Wt* ^ ■loCe«t*T.r./,r.T;))’*'’ (1.10) 


Ihe dependence of threshold (critical) drift velocity Uj, on the tenperature ratio 
Tg/Ti for hydrogen plasma is presented in Figure 1.5 [19] . It is apparent from the 
figure that in a fairly broad region of tenperature (0.1 < Tg/Ti< 8) , the EIC insta- 
bility is excited earlier (i.e., with smaller longitudinal points) , than the ion- 


sound. 
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Figure 1.5. Dependence of Critical 
Velocity of Electron on 
Tenperature Ratio Tg/Ti 
for Ion-Sound (I) and 
Ion-Cyclotron (II) Insta- 
bility 


f I , 

Figure 1.6 frcm [19] constructs the dependences 


and k„/kj^ on Tg/T^ for hydrogen plasma (approximately the same will be for plasma 
containing ions of any other type) . It follows from the Figure, for exanple, that 
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Figure 1.6. Critical Frequency Corresponding to 
the Boundary of Instability (^ 0 ) 0 ., 
Quantity = k2 and the 

Batio k^Aii Depending on Tg/T^ — 
Otenperature Ratios for Hydrogen 
Plasma. w/Wgj^ and correspond 
to tile same scale from the left, 
kj^Aji corresponds to the scale 
from the right. 


as Tg/Tj^ increases, the parallel ccxiponent of the wave vector k|| of the increasing 
waves rises (conponent along the field Bg) , In this case | ' . Figures 

1.5 and 1.6 refer to hydrogen plasma. Admixtures of heavier ions than change the 

threshold for developnent of ElC-instability. Thus, for example, in - 

plasma instability with develops at a lower threshold, even if the plasma 

is hydrogen for 80% of the ion content [19] . 


Vfe will sutrmarize the presented linear theory of ElC-instability. It is evi- 
dent tliat for studying the ElC-waves it is basically necessary to vise a kinetic 
equation, since the hydroi^namic theory of fluctuations in the magnetically-active 
plasma is insufficient, since it does not suitably take into consideration the 
thermal movement of the plasma particles. Because of consideration of thermal 
movement of plasma particles, a nvorber of new branches for oscillations in the 
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magnetically active plasma ^pear in the kinetic ^roximation 127, 26] ; weakly- 
and strongly-damping, long- and short-wave, and all of them are missing in the 
oold plasma. In particular, an electrostatic ion-cyclotron wave mode of Drurmond- 
Itosenblutli appears with frequency ajsnujg^ and with wave vectors directed at angles 
close to Tf/2 in relation to the external magnetic field. 

Ilieoretical stu<ty of the ElC-instability went beyond the framework of discus 
sion of a laboratory ‘'jcperiment (to a considerable degree because of the work of 
Kindel and Kennel [19]), and its results began to be used more often to interpret 
phenomena in the auroral zone. 

The develqpm^t of ElC-instability at altitudes of tlie upper ionosphere, for 
calm ccxiditions of the iono^ere, does not require very great density of the longi- 
tudinal currents (10® - 10^® electr./cm^ x s, or 2 x 10"^® - lO""® (19]) . Meas- 

urements in the auroral zone on satellites show the presence of longitudinal curraits 
vMch are quite sufficient for oscillations of the ElC-waves (see for exanple [47] ) . 

It is important to note that since the magnetic field is nonuniform because of 
longitudinal currents, the requirement of pressure balance results in thfs require- 
ment of pressure gradient. The pr^jssure gradient governs the drift instability with 
low threshold value of the current [30] . This effect is important with 
« 10“^ (As (i.e., for exaitple, in the region of diurnal polar cusp) , 

and consequently, is not suitable for the conditions of the upper ionosphere. 

The ElC-waves themselves have been observed in the auroral zone [24] , and 
simultaneously with ascending ions [44] . However, information on precisely vdiicli 
mechanism oscillates the ElC-waves (icn bundles or drifting electrons) has not been 
successfully obtained from the measurements [44] . If one considers the source of 
energy for the EIC waves to be ion bundles, then the question of obtaining these 
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bundles is still unresolved. 


Because of the lack of an unequivocal solution to the question of the source 
of free energy for ElC-wave oscillation, it is useful to also examine the possible 
mechanism for generation of ElC-instability not in the plasma with current, but in 
the presence of counter flows of ions [31] . 

In the plasma layer of the magnetosEiieric tail, the development of the explo- 
sive phase of a subsborm is aooonpanied by the appearance of hot plasma moving along 
the layer with average-mass velocity exceeding 600 - 1000 krt/s (32 , 33] . It is 
believed that these streams form in the region of fonnation of a neutral point at 
geocentric distances from the earth less than 10 Injection of plasma to the ^23 
earth among the geomagnetic field with velocities considerably exceeding the velo= 
city of the transverse drift, can cause fonnation of a oounterflcw of plasma as a 
consequence of reflection of the first stream from the nearer points in the iono- 
sphere. Under certain conditions, it is possible that different unstable modes 
will develop in this system, and in particular, electrostatic ion-cyclotron whicli 
is a probable cause of the development of anomalous resistivity in the auroral tube 
[15] . The results of satellite measurements [34] show that quasielectrostatic 
noise appecurs in the area of the plasma layer ( as well as in the auroral tubes) 
simultaneously with recording of longitudinal proton streams. 

Thus, the problem develops of determining the conditions for the development 
of ElC-w/ives in the plasma which is characterized by the presence along a uniform 
itagnetic field of two oounterstreams of hot '^anisotropic" protons (T^ ^ Tj^^) in 
the presence of "ariisotropic" electrons (T|| g ^ Tj^^) (witiiout consideration for 
the drift velocity of electrons along the magnetic field) . Since the distribution 
of energetic electrons over the angles within the plasma Ic^er remained isotropic 
at the moment the ejqjlosive phase of the substorm begins [48] one can therefore 
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) 

consider T^u “ 

Then the oc»idition ciitaijied in [31] for oscillation of the EIG-wave for the 
first harmonics looks like: 

■t 

' < (.-vVr.WM,* (1.11) 

v^iere V — drift velocity of the proton stream. Ihe inequality (1.11) has the following 
physical meaning. Increase in the stream velocity of the protons along the magnetic 
field diminishes the threshold for development of EIC-<lnstability. The parameters 
of tlie magnetospheric plasma during development of tiie explosive phase of the sub- 
storm satisfy., tlie inequality (1.11) [31, 36] . Thus, the observed ElC-waves can 
oscillate in the magnetospheric plasma during ?;nteractiai of the oounterstreams of /24 
hot protons with great perpendicular tenperature . 


Because of the inportance of taking into consideration the tenperature aniso- 
tropy in magnetospheric applications (for exanple, in the upper ionosEhere where 
the plasma is collision-free) , we will present formulas for this case [29] . Expres- 
sion (1.8) for tlie frequency of ElC-fluctuations is transformed into 


For the incretent of ElC-instability, instead of (1.9) we have 


( 1 . 12 ) 


J^N - 




''Jin « 






♦/a 


Respectively, instead of 



(1.10) we obtain 


T>h ’ 
Tlili *fli 



(1.13) 


(1.14) 
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1.3 , EVOimiON OP EIC-INSTABILITlf 


The general plan for evolution of ElC-instability excited by electron drift 
along the magnetic field can be presented as shown in Figure 1.7 [14] . 


p H-Electuxxv energy of drift I 

i 
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Ljs|^a ^rmal energy of electrons! ‘ r>f f»W-j 


^Nonoonerencj — jXon tliermai erterSyT 
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"Ion convection 
ses of ion energy! 


tron energy 


Figure 1.7. Schematic Illustratioi of Evolution of ElC-Instability [14] 

We note first of all that in its nature, the ElC-instability is weakly 
convective [21] . This circvmstanoe is important for inhomogeneous conditions of 
the ipper ionosphere where waves reach considerable anplitudes in the region of 
instability before wave convection from the instability iregion becomes important. 
In the regiai of maximum instability, the group velocity of waves transA'erse to 
the magnetic field approaches zero, 'idiile the grotp velocity along the magnetic 
field remains a finite quantity [21] . 

It is easy to show that by using (1.2) and the formula to determine 

■* ). ‘ slightly convective nature 

of the ElC-waves permits the instability to reach large anplitudes. In this case 
the qucisilinear effect of formation of a. "plateau" on the distribution function 
of electrons vMch results in relaxation of the oscillations, in the opinion of 
the authors [20, 22] cannot play its role because of the spatial effects and 




t 


t 


longitudinal electric fields occurring in the magnetosphere. Direct measurewents 
on the satellite S 3 - 3 [15] also indicate that ElC-inst^ility is saturated at 
a higher level than predicted by the quasilinear theory [1] . Ohe aitplitude of the 
ElC-^aves therefore rises until nonlinear mechanisms of saturation begin to cerate. 
One of these media \sms is broadening of the ion resonances [4] . In the opinion of 
the authors of [22] , in addition to broadening the ion resonances, the mechanism of 
longitudinal capture of electrons is active, vhen a large nuntoer of electrons are 
trajped, between the potential "walls'* of the wave. This, in turn, results in anoroar 
lous resistivity to the current of electrons. The theory of broadening of ion resai- 
ances is based on the concept that as instability develcps, plasma becomes turbulent. 
Orbits of ions moving in the turbulent plasma experience nonocherent perturbation 
because of the interaction of ions with stochastically changing caqponents of the 
wave field. Consideration for this interaction results in distortion of the condi- 
tion for particle resonance with the wave . This in turn, 

increases the niirtoer of particles which exchange energy with the wave. 

\\ 

Increase in wave damping on the ions results in replacement of the linear incre- 
ment hy the nonlinear This equals the linear increment supplemented by the 

term which corresponds to nonlinear damping Au , i.e., 


-*A.vo 


(1.15) 


The quantity Aai is associated with wave amplitude (5n/n) by the formula [4] 

/ ^ r f, A \ 


(1.16) 


vhere 




.Ik L { -=7 
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In this ceuMii the llniltat.^ion of wavs growth is determined idle condition 

I WW I 

t *0 » J (i,X7) 

The quantitative characteristics for restricticm in wave growth is the level of 
tujdbulence • ^l/nTt , where %-^-density of energy of electrical field of 
turbulent pulsations averaged for time. If in the expressicm for linear increment 
I jf^i ^^4 (1*9) , we ignore the term for resonance damping of a wave on ions, then 


the estimate of the turijulent level according to (1.16) and (1.17) is determined by 
the formula [4, 37, 12] ; . 


Vlt 






2 ci » 0.4 


(^y 


(1.18) 


Experimental confirmation of the nonlinear stabilization of ElC-instability according 
to the Dun-Dupree mechanism [4] was obtained in laboratory experiments (16, 18) . 


The stage of nonlinear limitation in the growth of ElC-instability is followed 
ky a stage of turbulent heating where the energy is tr^^<nsmitted fecm tl)e wave to v 
the particles. Wave energy density averaged for time (electrostatic plus Icinetic) 
for the first hannonics pouals fl] : 







(1.19) 


/28 


vhere , -Uc /lOp, , 1, —Debye radius of electrons. Tte first term in (1.19) 
corresponds to the electrostatic energy of the wave, the next two terms correspond 
to the ikinetic energy of electrons and ions respectively. Just as, usueilly in the 
magnetosphere kX© If so the kinetic energy of the particles is considerably 
greater than the energy of the field, while the ion kinetic energy is times 
greater than the kinetic energy of the electrons. Assuming that E] 5 ~ exp (yt) , 
the rate of energy transmission from directed movement of the electrons to kinetic 
energy of icxis participating in wave movement, is written as 


We will now eocamin® the process of energy transition from llC-imves to 
tlierroal energy of ions in the framework of the quaailinear theoj^ [9, 14] , We 
will define the heating rate as , 

^ vvi) ^ ^ (i‘2i) 

We initially focus attention on the "perpendicular” heating governed ty resonance 


particles, i.e. 






^ ^ <iV I h{Ni vf 


(1.22) 


where the equation for the distribution function of ions fj[ looks like; 

ai 7 s ***: ^ Vi, av* ) (v''” k. ♦ 


(1.23) 




Here %-frequency of node which corresponds to the linear theory; Jj^~Besselian /29 
function. Further, the ooniwtations use thd:^ Maxwellian distribotion function f^ 
and take into consideration only the first Imrmonics (n « 1) . 


By using (1.23), instead of (1.22) we obtain 


^ Jv. iCiXo-v.. 

9 % 2 Vi*- ** UV 

» r^Vx.V. ld*|4i, 


(1.24) 


Ey integrating for v, using 
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whare /^4 f *’i ^T; , 

E 5 ^ using th© ejqpression for the term of damping on ions in (1.9) , we obtain (for 
the average jc) ; 

T; 1 Ji " 

9^ * tjlc (1.26) 

Siitrlarly, for "parallel" heating 




(1.27) 


Where 0.15 for Ti~Te, u al.5. /3£ 

Hius, it is ^^arent from (1.26) and (1.27) that resonance heating of ions in a 
direction parallel to tlie magnetic field is an order smaller than the resonanoe heat- 
ing of ions in a perpendicular direction. 

We will now examine the nonresonance interaction of ions and waves for wliich 
instead of *hVi*) it is necessary to write the following in 

(1.23) 

vdiere y increment in mode correeponding to the linear theory . With integration for 
Yl in the region ^ (U»'‘ v£*^ |» »/ we obtain for turbulent heating of non- 

resonance ions in a perpendicular direction (for average K) : 
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%u ''•miEuy (1.29) 

Tuii)ulent heating of nonreaonanoe ions in a parallel direction (along tiie magnetic 
field) is exponentially small [14] . 

Xn this case one should recall that tiie effect of lieatj,ng nonresonanoe particles 
is ^parent, in the sense that (1.2B) describes the respcxise of tiie nonresonance 
particles to a change in wave anplitude. For exanple, v4ien the wave anplitude in- 
creases, the oscillatory kinetic energy of the particles associated with the wave 
also rises, and the nonresonanoe particles seem to bej heated [9] . 

Another situation develops vhsn ^e resonance conditions 
are distorted because of the interaction of particles with turbulent pulsations [17] , 
liien the 6- function is replaced fcy 



(1,30) 


> tn. ♦ ^ 

( u> - tfctr, - ♦ Cif * ♦ i* PO 


where — oonplete increment of wave increase which equals zero in tlie 

state of saturated instability; Dj_— -transverse diffusion term. In the region 

(i.e., in the resonance region) it is necessary to use the 
6- function. Far from the resonance region, where v)^ , it is 

necessary to take (1.30) , and is replaced by the average value [17] 

AlO, 

As a result, for ElC-instability in a state of saturation (with regard for (1.29)) 
one can obtain 

; (1-32) 


Formula (1.32) by magnitude yields the same rate of heating as the formula suggested 



However, (1.32) has a more axmon nature, since it is suitable for tlie case of 
wave anplitude saturation not only because of broadening of the ion resonances 
(formula (1,33) was obtained on tliis assumption) , but also for possible saturation 
of the wave aiiplitudes with the help of convection effects in limited plasma. 

Thus, one can state that distortion of the resonance conditions between veloc- 
ities of the particles and phase velocity of the waves makes it possible for the 
greater part of the ion distribution function to exchange energy with tlie electro- 
static ion-cyclotron fields. As a result, the rate of increase in ion tenperature 
is proportional to complete linear increment y^. If we only examined the resonance 

Tj 

particles, then the rate of heating ions' "j. “ would be proportional to |y ^ |. 

We will new dwell on the question of tlie mechanisms for saturation of turbu- /32 
lent heating for ElC-turbulence. Since the first ejqjeriments which studied turbu- 
lent heating and its consequences were made under laboratory conditions with fairly 
small dimensions of the xinit, therefore the corresponding theoretical constructions 
adopted the convective model for the heat emission mechanism, i.e., for saturation 
of turbulent heating. Correspondingly, the equation of thermal balance Icoked like 



L — dirension of the system. A similar approach in applying turbulent heating for 
saturation in longitudinal magnetospheric currents was develofjed in [37] , vhere tlie 
following formula (1.33) was obtained for the velocity of turbulent heating of ions 

» 4 V** , • It was considered in this case that the ac±ual ElC-insta- 

Tt » 4 H / 

bility is saturated by the mechanism for broadening of the ion risonances [4] . The 
quantity diminishes with a rise in ion tenperature (see (1.9)) . Since insta- 
bility occurs when u>u^, see (1.14), vbile u<-.f correspondingly, rises With an increase 
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in Ti, oonsequently, with a rise in icai tenf)erature '►u, while »^-*0 • 
Ihus, the effect of turbulent heating limits itself. In the condition on ttie boun- 
dary of stability, the "residual” turbulence heats the ions with sudi a rate that 
the processes of heat emission are capable of removing energj/ from the region of 
turbulence. In this case, one can write the equation for energy balance [37] : 


ttadtexi-X 



It should be stated that in the absence of heat emission, since turbulent heating 
in (1.33) is very rapid, ElC-instability will result in flashes of turbulence and /33 

local "spots" of liot ions [37, 38]. Ihus ElC-instability under certain cases can 
play the role similar to the Bunemanovskiy instability, result in brief flashes of 
ancnralcus resistivity, and not lengthy quasistationary states. 

Hie mechanism for heat emission v\hich guarantees saturation of turbulent heating, 
may be different [13] , vAiere energy from juitps to through the process of ion- 
ion collisions. In this case, if we consider that Tejj = const, while the thresliold 
of instability Ug ^preaches u with an increase in because of "perpendicular" 

heating of ions (see (1.14)) , then from an equation of balance for this mecnanism 
we can obtain a level of turbulence for the second stage of evolution of ElC-insta- 
bility, i.e., for turiDulent "perpendicular heating [13]: 

/ E. - ^ ■’‘u ■*<■>./ a-34) 

TJf. 

In this case ^ 5* • Ei^riments on S 3 - 3 [24, 15] forced us once again to 
attentively analyze the question of mechanisms for saturation of turbulent heating. 

In this case consideration of anisotre^y of the teirperatures of ions and electrons 
was basically inportant. The level of turbulence at the first stage of evolution, 
as shown above, is determined by formula (1.18) ; 
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licwever, [15] drew the conclusion that the second stage of evolution ("perpendic- 
ular” turbulent heating) does not significantly diange the quantity This conclu- 
sion was drawn based on measuremants on the satellite S 3 - 3, whidi showed that the 
measured level of turbulence e is close to ei, while was confuted from measure- 
ments on the same S 3 - 3 aocording to formula (1.18) . Therefore, decrease in the 
level of ElC-turbulence in the beginning of tiie second stage, viiich was expected 
aocording to the idea of "residual" turbulence [37] and according to formula (1.34) 
does not occur. Uius it was assumed in [15] that ’ The physical inter- 

pretaticMi of the conclusion that e ^^2 is based on the fact that [15] the ratio 
Tgjj during turbulent heating remains almost constant, since both teirperatures rise 
approximately at the same rate, so that u^ [see (1.14)] approaches u very slowly. 

Correspondingly, tlie level of turbulence ^2 slowly ^preaches a certain value 

00 

^2 vrfiich it adopts on the boundary of stability. Moreover, if in fontiula (1.34) the 
quantity vii is replaced by , where [11] f tlien one can obtain the 

quantity r-2 on the order of which coincides with tlie measured level of turiaulence 
[15]. Thus, tlie authors of [15] draw the conclusion that decrease in the level ^2 
as compared to is not a sufficient and significant circumstance for establishing 
the fineil level of turbulence. 

In concluding the section on evolution of ElC-instability , one can make the 
following sunmarizations . ElC-instability evolution is customarily divided into 
several stages. The first stage is frean the development of instability to its non- 
linear saturation at the corresponding level £ 2 ^. Laboratory ej^riments vhich have 
studied ElC-instability excited both ty currents [16, 18] , and by ion bundles [23] 
have shown that saturation of wave amplitudes occurs according to the mechanism of 
broadening of the ion resonances [4] . Then follows the seoond stage of evolution 
associated with t\nbalent h^iting. At tliis stage "perpendicular" heating of ions 
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occurs. Corresponding laboratory neasuronents [17, 14, 38] prove the existence 

of "perpendicular" heating of ions and agreement with theory. Because of heating 

of ions, tlie threshold of instability u,, increases almost to u, and consequently, 

the system approaches the boundary of stability (the condition of instability u> u , 

o 

whexre u is the drift velocity of electrons along the magnetic field) . This self- 
restriction of instability where the oonplete linear increment becomes very small 
with a lower level of turbulence G 2 , correspondingly results in a slowing down of 
the rate of turbulent heating. We note, however, that experiments on the satellite 
S 3 - 3 indicate the possibility of such a situation vbere E 2 . The final stat- 
ionary turbulent state is further defined from the energy balance in which the velo- 
city of turbulent heating in a state close to the boundary of stability must be 
balanced by the effective rate of heatlemission. i^ropriate laboratory experiments 
[18] demonstrate considerable transfer of energy and particles perpendicularly to 
the magnetic field, so that around the region with tiie current (i.e., around the 
filament) a dense ring of hot ions is formed. In the actual filament, the density 
in the nuittoer of particles is lower than in the surrounding ring of ions. The temper- 
ature of ions within the filament is such that the gyroradius of ions corresponding 
to it is on the order of or greater than cross dimensions of the fila- 

ment. Ihese same ejqjeriments indicate the appearance of suprathemal electrons dur- 
ing instability in agreement with the measurements made on the satellite S 3 - 3 
[15] , and the presence of energetic electrons strongly depends on the magnetic 
field in accordance with the theory of ancmalous resistivity developed in [11] . 

The next section covers an examination of approaches to tlie prcblem of ancrra- 
lous resistivity governed by ElC-tinbulence. 
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Since the reason for ElC-instability is exceeding of a certain critical value 
of current, then one can be convinced from general considerations that evolution 
of instability must result in a mechanism vdiich limits growth of the current. In 
other words, an additional electrical resistivity must appear which is called anoma- 
lous. 


Ihe standard formula for plasma resistivity contains the frequency 

of collision (v) of electrons with the scattering centers (ions, neutral atoms) in 
relation to pulse loss. If the plasma electrons build up certain types of oscilla- 
ticms for waves as a consequence of instability, then anomalous loss of pulse oocnirs 
(its transirassion to oscillations, i.e., collective movements of ions) . This loss 
of pulse is characterized no longer by v , but by a certain frequency v * . We will 
determine v* according to the formula [14] : 


^ v„ ^ (1 

vhere for the veLocity of change in electron distribution function fg because of 
interaction with waves one can write the following in the framework of the quasi- 
linear theory [9] : 



For ElC-waves <•< , and by leaving in (1.35) only the term with n = 0, 

after integration for Vj_ we obtain 
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wiiere Z(x) —tabular function associated with the Kranp function detenninsd above 
(W^x) by the correlation j H (k) t W(x'5 . ^ further consi^r that 

’ ^ for i . ; that f-jr- 1 . jf^ « , and by exam- 

ining only the first ion-cyclotron harmonics in oj|^, we obtained (using the asymptotics 
Z for small arguments) : 


hit-hf U 
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The first term in the brackets provides a contribution from the resonance particles. 

It is an order greater than the contribution frcxn the nonresonanoe particles v^iich 

0 




describes the remaining terms in the brackets, s,%:e 0.15 for T^"* Ti (14] . 

By leaving only the first term in tlie brackets, we obtain 




<ti> 

5i 

Finally, for the average Ewe have: 


vhere e --level of turbulence, vhile corresponds to the first term in (1.9) . 


Consequently, the effective quantity of anomalous resistivity [14] ; 

C hIT 


affl M — ■ — hi 


We note that turbulent "perpendicular” heating of ions mainly occurs because of 


nc;nresonance particles, «<(hile the effective frequency of collisions v* vhich deter- 
mines the pulse loss during collective collisions is governed by resonance particles. 
Ihe quantity of anomalous resistivity for ElC-instability which is in a state of 
saturation (according to the mechanism for broadening of ion resonances [4] ) with 
level 0.1 (wg£/a)p£)2 was conputed in [12]. Ihe authors started from the first 
racment in the equation 





r. *D* 


vhich yields 



(1.36) 


Here Df — diffusion coefficient in the spatse of velocities, e 3 _-''mentioned level of 
turbulence. Since (see for example [12f) ; 
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viiere — ^Besselian function on the order of n, v^ile 


A 

lO - r„-vr,, - * 





tO - Kntf,, - 
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then, in order not to solve the integral equation for D|, tlie resonance function of 
electrons was approximated by the Gaussian 

a-i*' 

where ^x/2 — half-width of resonance, Vc”critical velocity at which the diffusion 
coefficient has the maximum i i Co * ^^.Vc ) *• (. C''* 

After placing an expression for D® in (1.36) and integrating for velocities on the 


assumption of Maxwellian distribution for fg, in [12] Uie following was obtained 
for the quantily of anomalous resistivity 


and consequently. 




i. 




The same quantity v* in a quasilinear ^proximation was obtained in an earlier work 
[49]. 


The ion-cyclotron mode of Druimond-Rosenbluth as a .claimant for the main role 
in tlie problem of ancrralous resistivity of the magnetosphere obtained pcwerful si:p- 
port after publication of the measurement results on the satellite S 3 - 3 [15, 24] . 
These measurements proved the presence of the attained high level of ElC-turbulence 

• t// kTnT« * ^ 
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Oliis permitted the authors of [15] to assert that in the observed regions of EIC- 

turbulenoe, the "plateau" on the distribution function of electrons is not formed, 

«» • 

since for the measured level, the ratio . ' , vrtiile the quasilinear tlieory 

predicts Anomalous resistivity based on measured parameters of turbu- 

lence was ccjnputed in [15] , based on the law for pulse preservation in the electron- 
wave ^stem 


where Gq— imaginary part of electron dielectric function [28] . Adopting an ^rox- 

imation for random ftese during integration and monochromatic spectrum with k = E, 

• • I 

one can obtain •i* * . xC.I.t, E| / IHWtwu / where integration for d^Bd was 
done from ^ ^ using the value (see (1.2)) , we have: 


'!>♦* srr: 






( 


After taking, according to the measurements. 


’ JS- „ 



k./«i -o.i , tc. ^ S' -lo'' 
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Ji. . 

= 10"^ u = 6 X 10' cto/s, we find that • , while corresponcfling 

ancmralous resistivity n si x 10 ^ (Sim x m (vAiile the resistivity governed by part- 
icle collision ~ 10”3 ohm x m) . 


Field intensity governed by ancn^lous resistivity E,j~ 2 x 10"^ V/m. It should 
be said that these estimates in [15] were made on the assumption that the ion-sound 
turbulence is missing, altliough with u ~v^, the ion-sound or Bunemahovskiy turbulence 
must be pertuibed. Ibis circumstance was infli'-lt^sd in publications [39, 40] . In 
addition, accuracy of the estimate of the cof vhich is included in the 

formulas is not great and the result can easily toe exaggerated or underestimated by 
art order. 
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The obtai n ed fonrulas for anoTfitlous resistivity viere based on mechanisms of 
saturation vAiich did not take into consideration the boundaries of the field tubes 
in the magnetosphere. Since, according to current opinion (see for cxanple [40]), 
microwa<P structure of the longitudinal currents is possible, we will present an esti- 
mate of anomalous resistivity from [37] vAiich takes into consideration the boundaries. 

Prom the equation ^ vAiere t* '^*‘*‘*/*'*^ , 

vhile L~^imensiOn of field tube, the 

folk».ing was obtained ' • Mcmalous resistivity in this case 

according to the estimates in [37] is only an order greater than the quantity of 
resistivi^ because of particle collisions. 


2. MODEL OP m'LItffiD-Cycmi™ /4^ 

DajBLE POTENTIAL LAYER 

2.1. CWE-DlMENSlCJNAL MODEL 

OT inojINed layer 


We will examine the following problem. Assume that the electric current flews 
upwards along a uniform magnetic field directed along tlie axis Z (Figure 2.1) . We 
will consider that in the space Vq (i.e., above and to the right of the double layer) 
the directed velocity of thermal ions is negligible 


vhile, the electrons move downwards with fairly great directed velocity, ihe elec- 
trostatic potential on tlie vpper boundary of the layer equals <S>q. In the lower 
half-space (V^; , on the contrary, it is assumed that the electrons are fixed: 


( 2 . 2 ) 


vhile the ions are moving upwards with relatively great velocity. The potential on 
the lower boundary is assumed to be equal to* i. 
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Figure 2.1. Calculation of Inclined Double Layer 


CorrespOTdirigly, the indeK (gero) will fn^k the particles entering the leyr’er 
from above, vAiile the index (1) will indicate the particles entering it from below. 
It is assumed in this case that the parameters of the plasma and the field 
along the layer (i.e., along the y-axis in Figure 2.1) do not change ( d • o), 

Distribution of the electrostatic potential $ in the layer is determined by 
Poisson's equation: 

, V'tC'')* h.C?)] (2.3) 


(2.3a) 


If the layer is strong, then in the layer h;, « 1 

VtC?) vir« [«(, (x )i 

viiere n^^g — concentration of icai and electron, respectively equal to /44 

and fj^^g(u,r) — distribution functions of particles and velocities. Assuming that 
the distributicxi ption does not depend on the phase of the cyclotron rotation of 


I 


(2.4a) 


« 


the parhA-cle, we rewrite (2.4) in the fonn 

h-aiT 

where f (u,,/ Uj^)-- <listribution function of particla of corresponding type (the sispi 
i or e is emitted for brevity of writing) , satisfying the Vlasov equation, in a 
stationary case liaving the foliating appearenoe? 

.0 (2.5) 

As is known, solution to ei^tion (2.5) can be written in tlie form of a ranoom func-^ 
tion from integrals of particle novement (56] . Biis function is completely deter- 
minecl by assi,gning tlie boundary conditions. In order to search for the necessary 
movement integral in this case, (at least two) we will use tlie circimstanoe tiiat 
in the case where tlie transverse dimensions of the layer are much greater than the 


Larmor radius (we will be convinced of the correctness of this hypothesis later) , 
particle movement can be described in the framework of drift approximation. Ihen, 
as Swift showed [52), one can introduce the followi.ng movement integrals; 


and with consideration for the inertia drift of particles, 


(2.6a) 



(2.6b) /45 


■ u H 

where uh-.rinitial velocity of particle movement in the layer (u = Uq for electrons, 

Ui for ions), respectively <l>^^ = for electrons and $^=^l for ions, iDy- 

gyrofrequency of particles; " anfl ^ ^ 

tion of the leading center of tlte particle. By using the movement integrals (2.6a) 
and (2.6b) , we will write a solution to the equation (2.5) in the form 
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fact, by having the movement integrals as an argument, this functicm satisfies 
©q^tion (2,S) , and coinciding with she initial disuiibution function of frequency 
for velocity assigned on the layer bcwndaries, it also satisfies tlie boundary condi- 
tions. 

By substituting (2.7) in (2.4a), we detain ^ 

h »Z1T 5 p(^Uu , Ml) y* dUn (2.8) 

where the quantities Uj, and u„ are linked by the quantities uj? and u*j[ by the corre- 
lations (2.6) , from which, it follcws in particular that 


UidUi + 

uum 

(2.9a) 

au 

u,"aui: 




(2.9b) 


If the velocity of electric drift is small as compared to the velocity of longitud- 
inal movement u,i, expression (2.9b) is simplified and adqpts the following appear- /£6 


anoe; 


du. 


U»*c(u. 


(it. 1 \ 

^ ^ C»m) ’ 


or. [(«:)% 


(2.10a) 

(2.10b) 


uj is the velocity which particles would have if the electrical field was strictly 
parallel to the magnetic field. By substituting expressions (2.9a), and (2.10a) 
into (2.8), we find the particle concentration: 

Mu 
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Figure 2.2. Distribution of Electric 
Potential Transverse to 
Inclined Double Le^er in 
Cold Plasma 


vmich is the density of particles corresponding to the electrical field Ej_ = 0. 


It follows from expression (2.12b) and (2.10b) that 


F(u,r)fl(u,r 

» fevj.<t>S[uSp(UH”)/(u,r)'3 ciu,r 


(2.13) 


and the final expression for ion density looks like 




(2.14) 


' / >> V 2. 

As for tile electrons, because of their relatively sitall mass, the corrections 
associated mtli their inertia drift in a nonuniform electric field (in ei<pression 
(2.14), the terms proportional toy^^$ and Vj_$) are negligible and Ht, • (2.15) 

By substituting the ei^ressions cbtained for particle concentration in the /47 

Poisson's equation (2.3) , we obt.ain 

(2.16) 
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Taking into (consideration that in the exandned case t 






and ly introducing the Mven dielectric permeability ^ ■ i^ c'*/va where 
e)^ression (2.16) can be reMritten in the fom; 

KoMever, the cxxcrdinates x and Z, introduced as is shewn in Figure 2.1, are 
not very convenient for analysis of the examined task, since the bcxuidaries of the 
examined layer in this system are not coordinate lines. In tliis respect, we will 

'vV 

introduce a new Vcuriable \* €*’*•>•- / where a = tg0 (see Figure 2.1) . In these 
coordinates the boundosries of the layer are r^resented by the planes n = 0 and n = 


In this case, the final expression for Poisson's ec^uation Icxck like; 

{ * . r. ' 4 'i. V • I 


or 




1*1 


Hr ) 


(2.18) 


(2.18a) 


vhere 




<U‘ j <h»e «- 


Equation (2.18a) coincicles witli the known equation vhich describes the poten-- 
tial distribution in the (ioiible layer in the absence of an exteimal magnetic field 
(see for exairple, in equation (3.6) from [52J) with only the difference that in the 
ieft side of (2.18a) there is the nultiplier (1 +a 2e^) w'nich is usually iruch /4 8 

greater than a unit. Ihe quantity n^* and n^ is (conputed, as we have seen, as if 
the conponent of the electrical field vMcdi is transverse in relation to tiie kinetic 
field were missing, tnat is, in the same way as in the case of a direct double layer. 

This circumstance makes it possible, without repeating the mathematiccal computations, 
to use the results vhloh are known from an analysis of the <3irect potential layers 
in applicaticm to liie inclined double layer. In particular, if we assume = const 


40 


\\ 


for tlie roughest estijnates, then it is easy to see fioin (2.18a) that the thickness 
of the inclined layer along the aocis n is equal (with not very small a ) ; 

. h' lllTra =T>. 

from vMch it directly follows that the tliickness of the layer in tiie direction Z 
equals 

Ihe quantity oonprises^ as is knowri [57J, several Debye lengths: Dq = aXj^. By 
substituting the quantity Dq in e>q}ressic»i (2.20 ) , we find tliat the thidoiess of the 
inclined double Ic^er D » ap£, i.e., oomprises severul Larmor. radii of ions. Tlius, 
the hypothesis made above that the thickness of the double layer is great as oonpared 
to the Larmor radii of ions is ^parently fulfilled. 

However^ the numerical estimates presented above for the paraineters of the 
inclined double layer are very cpproximate . The fact is that in these calculations 
we assumed that the quantity was constant, ViAiile it significantly changes in a 

direction perpendicular to the layer. In this respect, we will examine the para- 

//"" 

. :: U 

meters of tlie layer and conditioiis for its ^istence in more detail. 


We will turn to the simplified hydrodynamic model examined in survey [50] in 
vdiich tenperature of the bundles of fast ions and fast electrons is assumed to be 
negligible (Tg = = 0) . In this case the distribution of fiiJX particles by 

longitudinal velocity is described by the delta-function, and equality (2.4a) is 


/4i> 


reduced to , 

( 2 . 21 ) 

By substituting the quantity Fj^^g from (2.21) into equality (2.11) , we obtain 
the concentration of ions in the form 


-// 


( 2 . 22 ) 


wiiere, as previously and 


Vt 


«. 

ti 


Ihe concentration of electrons, whose velocity of inertia drift is negligible, 


evidently equals 


Vl, -- H.O ( ^ , 


v^ere 


(2.23) 

V„« |v»i#Uc, 

By substituting the equalif'ies (2.22) and (2.23) in Poisson's equation (2.3) we 
can obtain an es^ression vdiidi detemines tlie distribution of the potential $ in the 
layer. licwever, the equation tlius obtained is still very coiplicated and is diffi- 
cult to analyze. In this respect, Sv/ift [52] sinplified the task even more, after 
assuming that the electrostatic charges are mainly concentrated near tne boundaries 
of the layer, while inside of it, a condition is fulfilled that is only correct in 
that case where tlie layer thickness is much greater than the Debye radius, ana there- 
fore is inapplicable in relation to the ion^-sound layer examined in the previous 
section. 


After equating ej^ressions (2.22) and (2.23), we find: 

t4>) 




(2.24) 


or in the dimensionless form: 


- ij 


(2.24a) 


iJere ^ is changed .in units and in units e‘I|i^ and distandes in the u^iits 

: r '>A^ere Pj-v:-Larmor radius of the ion with energy e$]^. 
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Equation (2.24a) is numerically integrated vdth assigned values and 
and fixed boundary conditions $ (o) » 0 with different values of (hii/r^) . As a 
solution, the soluticai is selected in v^^ch with » 1 , » o, iiie results of 

this integration with W^q = 0.1 and 0.001 are presented in Figure 2.2. On 

the x-axis the quantity^y^««^uS 9 , is Pitted, i.e., the distance (in units selected 
above) in a direction perpendicular the layer. We note that tlie quantity of the 
ratio (nj^i/n^) in this case was ec(tual to 0.6/ and not 0.1, as follcws from ti)e 

l( 

condition of Langmuir [50] . Ohis indicates the ii^licit approximateness of the 

Ccdculatiais. Ihe thickness of the layer, as is j^pparent from the figure, was equal 

\ 

to the two units indicated above, i.e., about i.5lip|. liiUs agrees fairly well with 

'■ '' ' ' 

the rough estimate presented above for this Quantity. 

; if ,,/• ' ■' ; ' 

In all the calculations presented above the tiMrperature of tlie bundle of ions 
and the bundle of elect^'ons was assumed to be negligible. As shown by Swift [54] , 
consideration for the final tenpsrature of the particles in tiie bundle does not 
change the main conclusion regarding the existence of solution (2.24a) for the 
random value The effect of the final temperature of the bundle is mainly reduced 
to increase in the thickness of the layer. Ihis only inproves the oonditions of 
^plicability of the discussed model based on the assuirption tliat tiie thickness of 
the layer is much greater than the Larmor ‘ radius of energetic ions. Ihe effect 
of tenperature of energetic ions on the width of the inclined l^er can be illus- 
trated ty Figure 2.3 [52] . It shows tlie distribution of potential in the layer in 
a model ^which contains a bundle of cold electrons above the layer and ions of tv® 
types! cold and liot, below the layer. As is apparent from the figure, incre^e in 
tenperature of liot ions actually results in oonsiderable (all the way to r/a 10) 
increase in the thickness of the layer. f 
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Hie cxsnditions for existence of the inclined double l^er in tlie ^prks of Swift, 
essentially liave not been studied. Hie fact is that in the model for cold bundles 
and cold plasma described above, tnese conditions are essentially reduced to the 
Langntuir oondition: 

In the hot plasma, the conditions for existence of the inclined double layer are 
reduced to the inequalities [50] 

Hie process of formation of the inclined double layer and conditions for its 
develcpnent have also not been studied in rJetail at all. Nevertheless, certain quali- 
tative considerations do exist on this account. Hius, ^ift, in the alrea^ cited 
publication [52] examined the conditions for build-np of the cyclotron (n = u) in- 
stability (from v^iich, the inclined double layer nost likely develops) , and shewed 
that the indicated instability develops in that case \diere the relative velocity of 
the electrons exceeds the threshold value 

vdiere a — factor on the order of a unit, vdiich is belcw the threshold of develcpnent 
of ion-sound instability [19] . 

If the plasma contains ions of two types, cold and hot, then the thresliold 
value of velocity for the bundle of electrons is equal to [52] : 


\Wt) \Titl 


(2.28) 




\\ 


vd»re andf' ^respectively the temperature and concentration of hot ions, and /53 
i\j— 'Concentration of electrrais. It is e^)parent from ei^ression (2.28) that if the 
density of cold ions is low (so that njj^® n^) , the threshold of development of tne 
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/52 


Figure 2.3. Distribution of Potential 
Transverse to Inclined 
Double Layer in Case of Hot 
Plasma 


examined instability drastically increases. Ohus, the presence of relatively cold 
ions is a necessary condition for development of ion-cyclotron instability, and 
correspondingly, the appearance of inclined double layers. 

In the model examined above for tlie inclined double layer, as in the model of 
a direct double layer, it was assuned that botii above the layer and below it, the 
plasma is uniform. In the real magnetosphere, tliis is not so of course. In partic- 
ular, in the exosphere, plasma density changes approximately according to the expo- 
nential law with a characteristic scale of altitudes h = T/big. In this respect 
Swift [52] assixnes th^t the area of possible existence of the inclined layer stretches 

along the force line of the geanagnetic field to a distance on the order of H. 
Below and above this region the force lines of the magnetic field can be considered 
equivalent, ^ssible configuration of the equipotentials is presented in Figure 
2.4 [52] . This same figure shews the distribution [along the meridian) of the 
energy of electrons accelerated in the l^er, as well as intensity of fhe elec- 
trical field at the level of the iono^here. / / 
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Figure 2.4. Possible Configurj^tlons of Equivalent Elec- 
tric Field in External Ionosphere and Distri- 
butions of Energy of Spilling Electrons Corres - 
ponding to it and Horizontal Electric Field in. 
the Ionosphere 


Figure 2.4a corresponds to the sirrpXest situation described above for the 
single inclined layer wliich is limited id altitude. As is epparent fron tlie figure, 
the field E in the ionosphere reaches the maximum in the center of the layer. 

Ihe energy of the accelerated electrons also has a characteristic latitudinal 
profile with maximum near the "Icwer" boundary of the layer. In this case, the 
characteristic scale of the picture in the meridianal direction, as v;s have seen, 
is 1 - 2 Larmor radii of energetic ions, that is, 1 - 2 km at the level of the 
ionosphere. 

Figure 2.4b presents two layers adjacent to each other in vdiich the horizontal 
component of the electric field is directed to opposite sides. 

Figvure 2.4c corresponds to a nultiplet structure of double layers and the 
regions of precipitation 'of accelerated particles vhich correEpond to them. 


Finally, Figure 2.4d correspondB to the case where outside the exanined layer 
there is an electric field with layer not associated and directed strictly pexpen- 
dicijlarly in relation to the magnetic field, so that the force lines of the latter 
are equipotential. 

The limitedness of extension of the area for particle acceleration in altitude, 
a^ respectively, along the meridian results in an interesting effect vdiidi is acoes-^ 
sible to. the ground dsserver. ^ fact is, as we have already stated, in the case 
of ions of an inclined double layer which are fixed in relaticxi to the ground observer 

th^ must move in relation to the ions ip«rards and downwards, and in relation to the 

(y)' 

observer to the north or to the scmth with a velocity determined by the Langmuir 
condition (2.25) . Ihis condition requires tliat the layer move alo^)g the force line 
with a velocity on the order of . On the condition of r^ni « 1 

and = 3 X 10^ cnv's (Wg = 2.5 keV) *this yields u^ = 2 x 10^ cnv^s. In this case, 
the velocity of movement of the layer (and correspondingly, the arcs of the aurora 
borealis) aloig the meridian equals u^/a , and with a ^10^, this velocity is on the 
order of 20 n)/s. Ihis does not contradict tiie experimental data. 

Ihese are the parameters of the inclined double layers vdiich are suggested 
bcused on the theoretical examination. Under laboratory conditions, these layers 
have not been studied as yet, so for experimental verification of the obtained laws, 
we will turn to the data of direct observations for thesd layers in the magneto- 
spheric plasma. 

2.2. DOUBIE POmNTIAL IN THE MAGNETOSPHERE; 

TWODIMENSIONAL MODEL OF IHE LAYER 

Figure 2.5 (lower panel of the figure) presents typical behavior of three 
corponents of the quasistetionary electric field vAien the satellite S 3 - 3 moving 
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at an altitude about 2000 km intersects the zone of the aurorae boreales [58] . 

Fran the data presented in the figure it is evident that the external ionosphere 
at invariant latitudes 68® - 72®, extremely inhomogeneous electric fields are 
observed with intensity in a direction transverse to the magnetic field (Ej^) up to 
300 rnv/m and in a longitudinal (E,|) up to 30 iilV/m, In this case, judging from the 
flight time of the satellite through the r^ion of the strong electric field, the 
boundaries of this region along the satellite trajectory comprise no more than 


several kilometers. 



Figure 2.5. Electric Field and Observed 
Stream of Electrcms During 
Typical Passage of Satellite 
S 3 - 3 Through Region witli 
Paired Electric Sliodk Vi^ves. 
Below is a more detailed time 
scan in a scale corresponding 
to the width of the arc. 
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Such small dimensions of the region for the existence of the field, like its 
high intensity , make it possible to identi:^ this region witli tlie double electro- 
static layer. Ihe quantities indicated above for the longitudinal and transverse 
components of the electric field indicate that this layer is inclined in relation 
to the magnetic field, and apparently oorre^nd to the mcdel described above for 
the inclired ion-cyclotron layer. This hypothesis is confirmed by the circumstance 
that the examined area, as can be seen from the data presented in Figure «.5, was 
obtadned in a broader region of intensive electrostatic turbulence. Ohis agrees /57 
quite oofipletely with the theoretical model of the layer suggested by Swift. 

Detailed analysis of tiie extensive material obtained on satellite S 3 - 3, 
permitted the authors of [58, 59] to draw the following conclusions regarding the 
parameters of the double layers in the magnetosphere: 

1. The double layers were found in the entire range of altitudes of action of satel- 
lite S 3 - 3 from 1000 to 8000 km, and the probability of appearance of tne layer 
increases with a rise in altitude. 

2. me cojuensions of tlie double layer along the meridian are (in calculation for 
the 100 km level) frcm 200 m to 10 km. 

3. The intensity of the transverse component of the electrical field in the layer 
reaches a quantity on the order of several tens of volts per mhter, all the way to 
1 V/m. Similar results were^ obtained previously from tlie data of launchings of 
barium clouds and streams at altitudes of several thousands of kilometers [60, 61} . 

The longitudinal component of the field averages an order weaker. 

4. The total difference in potentials transverse to the layer is several (all the 
way to 10 or even more) kilovolts . 

5. Double layers are chserved in the region of relatively intensive longitudinal 
currents Cjs (1 - 10) 10“^® Vcm^) / suid are acoompanied by plasma turbulence with 
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anplituae of the field fluctuations to 50 irtV/in. 

It is easy to see that the observed characteristics of the double layers in 
the magnetosphere on the v*»ole are fairly close to the properties predicted by 
theory. 

Hie difference in potentials existing transverse to the layer in this case is 
quite sufficient to accelerate the auroral electrons to the observed energies. Iliis 
circumstance makes it possible to hypothesize that the double l^ers observed in the ,/M 
magnetosphere must be linked to bright discrete forme of aurorae boreales. And in 
fact, direct comparison of ^)ecific cases of observed double layers with the aurorae 
boreales convincingly confirm this hypothesis [58] . Study of the link between the 
double layers in the magnetos£^re and the aurorae boreales can be significantly 
eipanded and deepened, as jpplied to precipitations of electrons of the "inverted" 

V type. 

Precipitations of this type were reported for the first time by Prank and Ackerson 
[62, 63] . The essence of this phenoroaion is the following. Figure 2.6 [64] shows 
the distribution (along the rocket trajectory) of luminescence intensity, density 
of the stream of energy of the precipitating electrons and their energy (in the 
region of the spectral peak) . As is apparent fron the figure, with intersection of 
the arc of luminescence, tlie energy of the electrons increases, reaches a certain 
neximum (in this case over 10 keV) , and then again diminishes, outliniiig on the 
graph a figure vAiich in shape is similar to an "inverted" letter v. In order to 
explain these features in the distributicxi of intensity of the stream and energy 
of the precipitating particles, Garnett [65] suggested that precipitations of the 
examined type are associated with the existence in the ipper ionosphere of an 
electric field vdiose characteristic distribution in space is described by a system 
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Figure 2.6. Distribution of Intensity 
of Precipitating Electrons 
and Brightness of luminescence 
According to Data of Satellite 
and Ground Stations 

of equipotentials of the type presented in Figure 2,7. In fact* it is easy to see 
that in this case the particles vMch move along the central axis of the figure 
intersect tlie greatest nuntjer of equipotentials* and correspondingly* tiie roaxiimm 
acceleration results. Bie energy of the precipitating particles diminishes the 
farther from tiie axis of the figure. 

We note that son«diat earlier, an analogous model of spatial distribution of 
the electric field in tlie auroral ionosphere and tlte equipotentials describing it 
was suggested by Carlquist and Bostrom [66] . Ihe model they suggested was based on 
results of observing the movement of inhomogeneities of aurorae boreales rays in the 
auroral arcs [66, 67] . Ihese observations vAiich shewed that irbeatogeneities ob- 
served to the soutii and to the north of the central axis of the arc move along the 
arc to different sides with a velocity required for the existence of the electric 
field " 0.5 V/m (in conversion for the level of the ionosEbere) . Since these 
intensive fields were never observed directly in tine ionosphere* the authors [66] 







original pack w 


also suggested a model sinilair to that deleted in Figura 2,7. 



Figure 2.7. Distribution of Equipotentials 
of Electric Field in Model "V- 
Shock Wave." 


Direct observations of the electric field in the magnetosphere that viere nade 
on the satellite S 3 - 3 shew that the double layers are encountered roost often in 
pairs or even series, and they are such that the electric field (E^^) in the adjacent 
layers is directed to different sides [58, 68] . 

Tlieoretical substantiation of the indicated distribution of electrostatic poten- 
tial was given by Swift [52, 54] . Ihe model he suggested is a develc^ent of the 
model exanvined above for a planf*? inclined dovible layer. It consists of the follo- 
wing. Assume that outside the region a^x < a, the longitudinal electric field is 
missing. Within the indicated interval, distribution of the potential at level z = 

Zq and at level z = Z]^ is assigned (similarly to Figure 2.1) : 

» 4”i (.*) 
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On the central axis of the figure, because of its suggested axisynitetry /61 

^1 .0 / ( 2 . 29 a ) 

Distribution of the potential in the layer, i.e., in toe region |xl£a and 
Zq <z < Zj^, is provided Poisson's equation (2.16) . Sino<b in the two-diirensional 
case we have examined, solution to this equation is excessively oomplicatea, we will 
introduce the following simplifications. Since, as follows from experixtent, 
we assume that everywhere with the exception of the axis, the layer (X = 0) ^C^ 


Tf 8 > ^ 


Then equation (2.16) adopts the appearance 

a**- ^ 9t Tta Vh,’* /’ 


(2.30) 


Concentration of particles included in the right side of equation (2.30) is 
defined by the untoown distribution of potential, and by the initial energy WgQ of 
fast electrons above the layer and of fast ions below the layer, as vjell as 
the distribution along the axis x of their concentration ngQ and However, random 

assignment of the last two quantities is impossible, since the correlation between 
themi aetermines the size of the jutp /which is already assigned. In 

this respect, Swift suggests the following procedure for solving equation (2.30) . 
we select a certain quantity where $2 is slightly less than $ 2 ^( 0 ) , 

and we numerically integrate in equation (2.30) . Integration is done until the 
equality 't’l W)* is fulfilled at a certain point X 2 < a. Hcwever, in 

order for the density of the charge to renf ^ *??«.nite everywhere, we additionally 
require that This condition can be satisfied only by the 

appropriate selection of the size of the ratio of primary density of fast ions in /62 
the region z <.zi to the density of fast electrons in tlie region z^Zq at the point 


After selecting tlie quantity R(x 2 ) in the necessary manner, we 
again integrate equation (2.30) on the hypothesis of linear relationship R(x) in 
the interval 0 <.X 2 . Me then select a new value for the potential on the axis 

of the figure ^^(0) < <i> 2 ( 0 ) and equation (2.30) is integrated with the sams boundary 
conditions to point X 3 at which A new value R(x 3 ) is selected. 

Ihis process is repeated many times until a sufficient number of profiles of /p (x) 
is obtained viiich correspond to different values of $( 0 ) , i.e. , integrating the 
central axis at different altitudes. The results or this calculation are presented 
in Figure 2.8a [54] . Ihe isolines of vdiich are thus obtained are presented in 
Figure 2.8b. Tlie, lower part of this fig\are shews the distribution along the x-axis 
of energy of electrons accelerated in the layer. It is apparent from the figure 
that both the configuration of tlie isollnes ^ , and the distribution of energy of the 
precipitating particles obtained in the model qualitatively coincide with similar 
ctiaracteristics of precipitations of the " inverted ‘‘ v typo, Tlie analysis made by 
Swift for the possible existence of the v-typt. double layer also contains, although 
in an irrplicit form, the answer to the question of vHnat governs the characteristic 
sliape of the layer. In fact, as we have just seen, a definite distribution of the 
quantity R(<)* transverse to the layer is necessary for the existence 

of a layer of the indicated shape. Vfe assume that the appearance of a double layer 
is governed by the intrusion of energetic electrons into a uniform (i.e., n^j^ = const) 
ionosphere. Then both the actual fact of existence of a v-type layer, and its posi- 
tion in space, and ^parently also the dimensions are determined by the distribution 
UgQ(x) . Thus, depending on the distribution of intensity of the bundle of primary 
electrons, both a v-shaped layer of random width (we recall, that in the case of a 
imiform bundle in uniform plasma, a flat layer is formed with width on the order of 
several Larraor radii) , and a simple layer of a certain configuration can appear 
(see Figure 2.4) . 
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lims, the ej^eriiiEntal data examined above make it possible to assert that 
double potential layers are just as widespread a phenomenon as anomalous resistivity, 
and the difference in potentials existing in them is canmensurate with the drop in 
voltage in the region of anomalous resistivity and is Sufficient to accelerate part- 
icles to typical auroral energies. 

It is thus apparently difficult at this time to state which of the discussed 
phenomena can be viewed as tlie main mechanism for accelerating auroral electrons 
along the force lines of the magnetic field. It should be noted however, that judging 
from the fact that double layers in all cases vhere tliey Were observed, were sub- 
merged into strongly turbulent plasma, as well as from the fact that tneir transverse 
dimensions are 1-2 orders smaller than the typical wididi of tiie zone of precipi- 
tation of energetic electrons, one can assume that on the] wlx>le they play a secondary 
role in tlie energetics of auroral phenomena, being responsible more likely for the 
relatively narrow and short-lived, although clear discret’e forms of aurorae boreales. 
Thus, double layers in the magnetospheric plasma are apparently a certain secondary 
effect or partial, although they may be a extreme manifestation of electrostatic 
turbulence of plasma in the presence of strong longitudinal currents. 
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